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We also show by explicit calcnlation that the hard, soft and jet functions in SCET are 
regularization-scheme independent. 
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1 Introduction 


Higher-order calculations in QCD result in loop integrals that are often ultraviolet (UV) 
and/or infrared (IR) divergent. The standard method to deal with these singularities is 
dimensional regularization, where space-time is shifted from 4 to D = 4 — 2e dimensions. 
The UV and IR singularities then manifest themselves as poles 1/e^. 

There are several variants of dimensional regularization. The most common scheme 
is conventional dimensional regularization (cdr), where all vector bosons are treated as 
D-dimensional. Prom a conceptual point of view this is the simplest possibility and guar¬ 
antees a consistent treatment. However, cdr has some disadvantages. Apart from break¬ 
ing supersymmetry, it is also not directly compatible with the helicity method and other 
computational techniques that rely on 4 dimensions and, hence, leads to more tedious 
expressions in intermediate steps of a calculation. Therefore, it is often advantageous to 
use other schemes, such as the ’t Hooft-Veltman scheme (hv) [1], dimensional reduction 
(dred) [2] or the four-dimensional helicity scheme (fdh) [3]. 

The result for a physical quantity such as a cross section is of course finite and must 
not depend on the regularization scheme that has been used. However, in practise such 
a result is obtained as a sum of several contributions, which usually are separately diver¬ 
gent. Therefore, these partial results can depend on the regularization scheme. It is often 
advantageous to use regularization schemes that are adapted to the technique used for the 
computation of a particular contribution. In order to be able to consistently combine the 
various partial results it is then imperative to have full control over the scheme dependence. 

The key observation is that the scheme dependence is actually intimately linked to 
the structure of UV and IR singularities. The singularity structure in fdh and dred is 
best understood if the (quasi) 4-dimensional gluons g are split into D-dimensional gluons 
g and = 2e scalars g. Prom a conceptual point of view these so-called e-scalars g can be 
treated as independent fields with an initially arbitrary multiplicity N^. The identification 
= 2e is to be made only at the end of a calculation. The decomposition of g into g and 
g has to be made in dred as well as in fdh. This seems to be a disadvantage of these 
schemes. However, it is useful to gain insight and to derive the scheme dependence, and 
for practical purposes, such an explicit separation is often not required. 

The contributions of the e-scalars are UV and IR divergent, resulting in terms of the 
form {N^y/e^. It is precisely these terms that - after setting = 2e - induce the scheme 
dependence in partial results. Por a physical cross section the poles in e have to cancel, 
including poles of the form N^/e. This entails that the scheme dependence for a (finite) 
physical result can be at most 0{N^e^) and, hence, will vanish in the limit e ^ 0. At 
next-to-leading order (NLO) this has been explicitly demonstrated [4]. However, virtual 
corrections generally are UV and IR divergent and, therefore, scheme dependent. To find 
this scheme dependence the structure of UV and IR singularities has to be understood for 
a gauge theory with gluons and e-scalars. 

Regarding the UV singularities, the main point is that treating the e-scalars as inde¬ 
pendent fields induces additional couplings. The independence of these couplings and their 
UV renormalization was already required in the equivalence proof of dred and CDR [5-7] 


2 


and in explicit multi-loop calculations in dred [8-10]. It has to be stressed that also in 
FDH the couplings have to be treated as independent [4, 11], 

The development regarding the scheme dependence related to the IR divergent part 
beyond NLO is more recent. The structure of the IR singularities for massless gauge 
amplitudes has a remarkably simple form [12-15]. It can be expressed in terms of the cusp 
anomalous dimension 7 cusp and the anomalous dimensions of the quark and gluon, 7 ^ and 
7 g, respectively. These anomalous dimensions have been extracted from explicit results of 
form factors computed in CDR and are consistent with other processes. 

It seems natural to assume that this structure can be extended to other schemes by 
applying the split of g into g and g. This results in modified (i.e. scheme dependent) 
anomalous dimensions. At NLO, this leads to results that are consistent with the well- 
known scheme dependence of NLO amplitudes [16]. Based on this assumption, 7 cusp) 7 q and 
7 g have been extracted in the fdh scheme at NNLO [17, 18], by comparing the generalized 
IR structure to explicit results of two-loop amplitudes for the 7 * —>■ qq and H ^ gg form 
factors and the process qq ^ g^- Considering all these processes together yields an over¬ 
constrained system for the extraction of 7cusp) 'Iq and 7^ in the fdh scheme. The fact that 
there is a solution to this system suggests that fdh is a well defined scheme beyond NLO. 

The main results of this paper are the following: First, we will provide further evidence 
that with a proper definition fdh can be used for loop calculations beyond NLO. To this 
end we show that the anomalous dimensions 7 cusp) Jq and 7 ^ can be computed directly 
in soft-collinear effective theory (SCET) [19-27] by relating them to the jet- and soft 
functions. We repeat the original calculation of the quark-jet function [28] and gluon-jet 
function [29] in the fdh scheme and also determine the soft function in fdh. This gives 
us an independent determination of 7 cusp) Jq and 7 ^ in the fdh scheme and the results 
we find are in agreement with previous findings. Note that the fdh as we use it [4, 17] is 
slightly different from previous implementations [30]. 

Second, we extend the scheme dependence study to dred. While the anomalous 
dimensions in DRED are the same as in fdh we also need to consider amplitudes with 
external e-scalars. Determination of the IR structure of these amplitudes requires the 
knowledge of 7 ^, the anomalous dimension of the e-scalar g. We compute 7 ^ in SCET via 
the calculation of the 5 -jet and soft functions and give the generalization of the IR structure 
to amplitudes with external g. Furthermore, we verify that this result for 7 ^ is in agreement 
with the result extracted from an explicit computation of 77 —>• 55 at NNLO [31]. We thus 
obtain a complete understanding of the relations between NNLO amplitudes with gluons 
and massless quarks computed in cdr, hv, fdh, and dred. 

Finally, we gain insights into how the regularization-scheme dependence cancels for 
fully differential cross sections at NNLO. While a complete study of this issue is beyond 
the scope of this work, our calculations in SCET show that the jet- and soft functions 
are separately scheme independent. The same is true for the hard function. Hence, if the 
cross section is written as a convolution of hard-, soft-, and jet functions it is manifestly 
regularization-scheme independent. Recently there has been a lot of activity in performing 
fully differential NNLO calculations using the SCET framework. This development started 
with the computation of top-quark decay [32] and has then been extended to more generic 
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cases [33-35] . The results of our work show how to apply a particular regularization scheme 
for the calculation of either the hard-, soft- or jet function. For each of these building 
blocks separately, the most convenient regularization scheme can be used. This opens up 
possibilities for further technical advances. 

The paper is organized as follows: In Section 2 we briefly review the various regu¬ 
larization schemes and discuss how they affect the IR structure of scattering amplitudes. 
Section 3 is devoted to the computation of the anomalous dimensions that are required for 
the IR structure. These computations are done in SCET. An alternative determination 
of the anomalous dimension of the e-scalar is presented in Section 4, where we extract 7 ^ 
from the gluon form factor computed in DRED. In Section 5 we use these results to obtain 
explicit transition rules for two-loop amplitudes between HV and fdh, as well as between 
FDH and DRED. The transition rules are then checked with explicit examples. Our con¬ 
clusions including a discussion on the scheme independence of cross sections at NNLO are 
presented in Section 6 . Finally, we give some explicit results of the SCET computations in 
Appendix A and list the required anomalous dimensions and (3 functions in all schemes in 
Appendix B. 

2 Schemes and structure of IR singularities 
2.1 Regularization schemes 

Dimensional reduction has been shown to be mathematically consistent [36] and equivalent 
to dimensional regularization [ 6 , 7] on the level of IR finite Green functions. In the way 
we define it, the fdh scheme has the same properties. The consistent implementation 
of the considered regularization schemes requires the introduction of three vector spaces. 
Apart from the strictly 4-dimensional space (4S) with metric two infinite-dimensional 
spaces have to be introduced, the quasi 4-dimensional space Q4S [36-38] with metric 
satisfying g^^ = 4 and quasi D-dimensional space QDS with metric g^'^ satisfying = D. 
The structure Q4S D QDS D 4S is reflected in the properties of the various metric tensors: 
g^'^gup = g% and g^^’^gup = g^p- 

For a detailed discussion and a precise definition of the four considered regularization 
schemes (rs) we refer to Ref [4]. Here we only repeat the most important aspects to 
facilitate the following discussion. The various RS differ in the way “internal gluons” (part 
of a one-particle irreducible loop diagram or unresolved final state gluon) and “external 
gluons” (all remaining gluons) are treated. This is summarized in Table 1 taken from 
Ref. [4]. Since external gluons are treated as stricly 4-dimensional in fdh and HV these 
schemes are best adapted to be used in connection with the helicity method. 

The cleanest way to understand the scheme differences is to consistently apply the 
split of the (quasi) 4-dimensional gluon into a D-dimensional gluon and an e-scalar. This 
is done at the level of the Lagrangian writing the field of the 4-dimensional gluon field of 
FDH and DRED as -I- A^, where A^ and A^ are the D-dimensional gauge field 

and the e-scalar field, respectively [5]. We will denote the associated ’particles’ as g and 
g, respectively. The e-scalars have an initially independent multiplicity and the metric 

associated with g satisfies the orthogonality relation g^'^g^p = 0 and g^’^g^u = 
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CDR 

HV 

DRED 

FDH 

internal gluon 





external gluon 

5^" 

5^" 




Table 1. Treatment of internal and external gluons in the four different RS, i.e. prescription 
for which metric tensor is to be used in propagator numerators and polarization sums. 


Scheme differences have their origin in UV and IR divergent contributions due to these e- 
scalars. These contributions are of the form {Nf^Y/e^ and after setting —)• 2e result in the 
scheme differences. This connection to UV and IR singular terms allows for a completely 
systematic treatment of the RS dependence. 

Regarding UV renormalization, fdh and dred behave in the same way. The possible 
split of internal gluons into gauge fields and e-scalars implies that in principle five different 
couplings need to be distinguished (see in particular [7, 8, 11]): the gauge coupling as, the 
gqq coupling a^, and three different independent quartic ^-couplings a/^s,i with i = 1,2,3. 
In general, we write the perturbative expansion of a RS-dependent quantity V^®({a}) as 




«e\"' f Ci4€,l\^ /«4e,3V .^^rs 

A, 


dvr/ 


dvr 


dvr 


dvr 


mnklj • 


Accordingly, the (3 functions for as and ae in full generality are written as 


2 ^ 

^ dg? dvr dvr 


E 


Vdvr 




E>2 

„ _ Q<e / Cks 

2 /I .ir 471- E-/ \ 47J- 

S>2 


d ap 


dg^ dvr 


dvr 

dvr 


Q:4e,l 

dvr 

Q:4e,l 

dTT 


04€,2Y /a4e,3V osRS 


dvr 


dvr 


A 


mnklj ’ 


*^4e,2 \ neRS 


dvr 


dvr 


A 


mnklj 


( 2 . 1 ) 


(2.2a) 

(2.2b) 


with analogous expansions for the /3 functions for a^p^i. In the sums, S > 2 is an abbrevi¬ 
ation for m + n + k + l + j > 2. The later results of the present paper will show that the (3 
functions of the are not needed and that we do not need to distinguish between them; 
hence we will often denote them generically by a^^Y Note that in Eq. (2.2) all quantities 
are finite and the scheme dependence is 0{Np). Thus, after setting —)■ 2e and then 
e —)■ 0, the scheme dependence disappears and we refrain from using an RS label on the 
l.h.s. of Eq. (2.2). In particular we write as and a^ without an RS label. 

According to Table 1, in dred external gluons are (quasi) d-dimensional. The de¬ 
composition of these external gluons into g and g also allows to avoid all problems related 
to factorization theorems [39] in DRED regularized QCD. However, this split results in a 
larger number of ’independent’ diagrams. Applying the decomposition of g into g and g 
then implies that in DRED amplitudes with external e-scalars have to be considered. This 
is not the case in the other schemes. As this leads to additional complications, we will 

^We remark that in practice the couplings can often be identified; only the bare couplings and the 
associated renormalization constants and functions must be kept different. Section 5 will provide further 
discussion and examples. 
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first restrict our discussion of the scheme dependence to the schemes CDR, hv and fdh in 
Section 2.2. Then we will consider bred in a second step in Section 2.3. 

2.2 IR structure in CDR, HV and FDH 

After UV renormalization, on-shell scattering amplitudes in massless QCD still contain 
IR poles 1/e^. In the framework of CDR it has been shown that these singularities can 
be subtracted in the MS scheme, using the procedure described in [12-15, 40-42], via a 
multiplicative renormalization factor Z which is a matrix in colour space. This can be 
generalized not only to the hv but also to the fdh and dred schemes [17, 18]. 

For the following discussion we find it more convenient to work with amplitudes 
squared. More precisely, we consider 

Vl«^*(e,V,,{p}) = 2ReMr(TA^oM)l^'^"*(e,iV.,M)), (2.3) 

where |.4^®*(e, {p})) is a UV renormalized, on-shell n-parton scattering amplitude con¬ 

taining IR poles and (.AQ®*(e, Vg, {p})| is the corresponding tree-level amplitude^ Both the 
e- and the V^-dependence differ in the four regularization schemes. For the moment we 
restrict ourselves to CDR, hv, fdh, as indicated by the label RS*. Then the regularized ex¬ 
ternal gluons behave completely as gauge fields and do not have to be split into gauge fields 
and e-scalars. The set {p} denotes the set of partons of the process under consideration 
and contains only quarks or gluons. 

The regularization-scheme dependence of is related to the IR poles and can 

be absorbed by a scheme-dependent factor We can define IR subtracted finite 

squared amplitudes as 

=2Re(^r(TiVoW)l(Z^"*(e,W,{p},/i))-'|.4^®*(e,W,{p})), 

(2.4) 

where p represents the factorization scale. The expression on the l.h.s of Eq. (2.4), 
denotes the hnite remainder of the amplitude where the poles have been subtracted in a 
minimal way. still depends on e (and W) but does not contain poles 1/e^ any longer. 

Hence, the limit e —0 can be taken and then we obtain a scheme independent hnite matrix 
element squared 


= lim A4fu®b(e, Ae,{p},/x). (2.5) 

(iV)e-i>0 

The limit (A)e —>■ 0 indicates that hrst we set —)• 2e and then e —)• 0. To put it differently, 

after setting —)■ 2e, the scheme dependence of is only in the terms 0{e). 

The starting point for a typical NNLO calculation is the computation of the two-loop 
virtual corrections in a particular regularization scheme. This corresponds to as 

dehned in Eq. (2.3). To understand the IR divergence structure and obtain transition rules 
between schemes we want to exploit the relation of the scheme-dependent Ad®"®* to the 

^Strictly speaking, the tree-level amplitudes in the RS*-schemes do not depend on Nevertheless, we 
keep the dependence on V in the notation to simplify the generalization to dred in Section 2.3. 
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scheme-independent A^gn- The key quantity for this is the scheme dependent factor 
to which we turn now. 

The all-order amplitude {p})) in Eq. (2.4) is independent of the factoriza¬ 

tion scale /i. It follows that the IR subtracted amplitude squared satisfies a renormalization 
group equation (RGE) 

{p}, p) = {p}, p) iV„ {p}, p ), (2.6) 

where the anomalous dimension r^®*(A^£, {p},p) is related to the factor through 

T^^*{N,,{p},p) = -{Z^^*{e,N,,{p},p))-^^^Z^^*{e,N,,{p},p). (2.7) 

This equation can be formally solved to obtain a path-ordered exponential with respect to 
colour matrices 

Z^^*{e,N„{p},p)=Vexp H ^T^^*{N„{p}, p'). (2.8) 

J 

In [12-15] it has been shown that in cdr the general structure of the anomalous 
dimension operator F, which controls the IR divergences of QCD scattering amplitudes, is 
exactly known up to two-loop level and only involves colour dipoles. In those papers it was 
also conjectured, by using soft-collinear factorization constraints and symmetry arguments, 
that this simple structure is more general and it is valid to all orders in perturbation theory. 
Generalizing this from CDR to other schemes and suppressing the dependence on N ^, we 
write according to Refs. [17, 18] 

r“-({p},M) = E *?P In + E 7.”-. (2.9) 

«,J) « i-1 

where Sij = ±2pj • pj + iO, the sign “-I-” is chosen when both momenta pi and pj are 
incoming or outgoing and the sign ” when one momentum is incoming and the other 
one outgoing. The first sum in Eq. (2.9) runs over all pairs i ^ j of distinct parton indices 
i,j € {1, 2,... , n}, where n is the number of external partons. The universal quantity 7 )[(fsp 
that appears as coefficient of the two-particle correlation term, Tj ■ Tj = T?T^, is called 
“cusp” anomalous dimension. The quantity yf®* is a single-particle term which depends 
on the type of the external particle, 7 ^®* = 7 ?'®* in the case of a (anti)quark and 7 ^®* in 
the case of a gluon. The explicit form of the colour generator associated to the f-th parton, 
T“, is as follows; Eor final-state quarks or initial-state antiquarks, the colour matrices T 
are defined by (T'^)^^^ = where is a SU(A^) generator. For final-state antiquarks or 
initial state quarks one has instead (T^)^^^ = —while for gluons (T'^)^^ = 

As a consequence the IR structure can be described by a set of three constants, which 
depend on the scheme 


RS* € {cdr, hv, fdh} : 


^ RS* _RS* _RS* 
TcuspjTg ) Tg 


( 2 . 10 ) 
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Thanks to the simple structure of the anomalous dimension matrix F, one can find an 
explicit solution for the perturbative expansion of Z. It is also possible to drop the path¬ 
ordering symbol in Eq. (2.8) since the colour structure of F is independent of /i. The 
following notation is often introduced 

5 ;a, ( 2 . 11 ) 

n ^ i 

where the last equality follows from colour conservation, Ci = Cq = Cq = Cp for 
(anti)quarks and Ci = Cg = Ca for gluons. 

All scheme-dependent quantities introduced so far potentially depend on all couplings 
{a(/i)} = {q!s(//), ae(/r), (//)}. Thus, in general the perturbative expansion is of the 

form of Eq. (2.1). 

Solving the differential equation Eq. (2.7) one obtains a perturbative expression for 
hiZ^®* which also depends on the /3 functions. Suppressing the arguments, in particular 
the dependence on the process {p}, it can be written up to NNLO as 


In Z^®* = 


a 

dvr 


■p/RS* ■pRS* 


4g2 


2e 



QORS* 

op 


mnklj 


■p/ RS* 

^ 1 


16e3 


^RS* ^ ■pRS* 
^mnklj * 1 

4? 


+ 


p/ RS* 
mnklj 

IGe^ 


+ 


pRS* 

mnklj 

4e 


+ 0{a^). 


( 2 . 12 ) 


Here the sum S = 2 denotes a sum over all terms satisfying m + n + k + l+ j = 2, and the 
following vector notation for terms involving pure one-loop quantities has been used: 


^ -pRS* _ -pRS* ■ 

a • 11 — ccs 1 j^oooo w ae . 


uioou 








^RS* ^ -p 
^mnklj ' 


RS* — 

1 — 


/QSRS* "pRS* _i_ /OeRS* "pRS* 

Pmnklj ^ 10000 Pmnklj ^ 01000 


I /o4e,lRS* -pRS* I n4e,2RS* -pRS* i n4e,3RS* -pRS* 

' Pmnklj 00100 ‘ Pmnklj 00010 ' Pmnklj 00001 ’ 


(2.13b) 


and analogously for the combinations involving T'^. The dependence of F on the individual 
couplings and the appearance of the different functions constitutes an important differ¬ 
ence to the CDR case, where only the a* and /3® terms appear. It can be obtained by setting 
cke, a 4 e,i 0 in Eq. (2.12) and identifying Fmoooo = etc. 

Eq. (2.12) shows that the one-loop IR divergences are described by the one-loop co¬ 
efficients of F', which depend on the process-independent quantity 7 )fus*p, and of F. Both 
anomalous dimensions depend on the partons involved in the process. At the two-loop 
level, the full 1/e^ and parts of the 1/e^ divergences are predicted by one-loop /3 and 
T coefficients. The remaining 1/e^ and the 1/e poles are described by genuine two-loop 
anomalous dimensions. 

Eq. (2.4) together with Eq. (2.12) allows to describe the RS dependence of the squared 
amplitude 
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• CDR-HV: Since internal gluons are treated in the same way in CDR and HV we have 

2 CDR _ 2 hv anomalous dimensions are the same in these two schemes. 

The difference in the squared matrix element comes entirely from using different 
metric tensors for the polarization sum due to external gluons. In cdr, where external 
gluons are D-dimensional, this polarization sum involves whereas in HV is 
to be used. 

• HV-FDH: Since internal gluons are treated differently in HV and fdh we have 
ZHv ^ z^dh ggg^ ^j^g anomalous dimensions are not the same in these two schemes. 
This results in further scheme differences of the squared matrix element. However, 
external gluons are treated in the same way in HV and fdh and the metric tensors 
in polarization sums are the same in the two schemes. 

2.3 IR structure in DRED 

Understanding the IR structure of dred processes with external gluons is more compli¬ 
cated. Each external quasi-4-dimensional gluon can be split into a g and a g, and the 
squared matrix element for a process with ^g external gluons can be decomposed into 2 ^^ 
terms. Following Ref. [4], we can write for the amplitude squared for such a process 

M°^^^{...gi...g#g...)= Y. ••• E (2.14) 

9i&{g,g} 5#9G{9.9} 

Reinstating all variables explicitly, we write the same relation in a more compact way as 

=^M°^^°(e,iV„{p},/r). (2.15) 

{p} 

Hence, the partons appearing in the list {p} on the r.h.s. can be either quarks or g, g, 
but not full quasi-4-dimensional gluons. We stress that practical calculations are not as 
complicated as implied by Eqs. (2.14) and (2.15). The l.h.s. will typically be computed 
directly as a whole with quasi 4-dimensional gluons, i.e. 4-dimensional numerator algebra. 
Even the renormalized couplings ag, 04 ^ can be identified, see section 5 for further 
discussion. However, from a conceptual point of view each term in the sum on the r.h.s. 
of Eqs. (2.14) and (2.15) can be considered as an independent process and the couplings 
as independent. Then, each of these processes behaves as the processes in CDR, HV, fdh 
discussed in the previous subsection, and it becomes possible to understand the IR structure 
and construct IR subtraction terms and transition rules to other schemes. 

For each process on the r.h.s. of Eqs. (2.14) and (2.15) a corresponding factor 
Z(e, Nf;, {p}, p) and a subtracted squared amplitude 7VI°,(if°(e, W, {i?}) A^) can be con¬ 
structed, like for in Eq. (2.3) and Eq. (2.4). Overall, one can then define the full 

subtracted squared amplitude in dred as 

M, m) =Y {p}, p) ■ (2.16) 

{p} 
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(2.17) 


It satisfies an equation analogous to Eq. (2.6), 

{p}, = /^)-^™(e, {p}, p ), 

{p} 

The for the individual parton sets {p} satisfy relations analogous to Eqs. (2.7), 

(2.8) and (2.9). Likewise, the subtraction factors Z can be written as 


InzDRED 







(2.18) 


O/QDRED 


■p/DRED 

1 


16e3 


ODRED . -pDRED 
^mnklj * 1 

4^2 


■p/DRED 

mnklj 

16e2 


■pDRED \ 


+ 0{a^). 


Like in the corresponding Eq. (2.12) the arguments are suppressed. An important difference 
to the RS* schemes is that in dred the individual split processes {p} have to be used. This 
implies that the set of 7 ’s needed to describe the IR structure is different in dred compared 
to the other schemes. 


DRED DRED DRED DRED 
'CUsp ’ iq ’ ig ’ Ig 


(2.19) 


This should be compared with Eq. (2.10). There are however several obvious relations, 
since internal gluons are treated equally in fdh and dred: 

( 2 . 20 a) 
( 2 . 20 b) 
( 2 . 20 c) 


- _ ^FDH _ DRED 

Tcusp — /cusp /cusp ) 


- _ „,FDH „DRED 

= 7, 


79 = 7^°'' = 7, 


Q 

DRED 
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Thus, the e-scalar anomalous dimension ^P^-ed jg only additional ingredient in dred. 
To highlight this, we introduce the notation % for this quantity, 


7c = 7r"° • (2-21) 

It is instructive to compare the individual processes with external g or g in dred to 
a process in fdh. The squared amplitude for a process with at least one external g has 
an overall factor iVc from the e-scalar polarization sum. As long as we consider the UV 
renormalized, but not yet IR subtracted matrix element, we cannot set (A^)e ^ 0 since 
there are still IR poles present. However, once these have been subtracted, the squared 
matrix element is free of poles in e and still contains a factor Hence, 


and 


...) = , lim ..~g...)=0 

{N)e^0 


( 2 . 22 ) 


lim 

{N),^0 


M^^^°{...gr...g#g...) 


, lim • • 51 • • • 5#9 ■ ■ •) = >lfin(- • • 5 • • •), 

(Ar)e-fO 

(2.23) 
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i.e. once the amplitudes are properly subtracted and the limit (N)^ —0 is taken, processes 
with external g do not contribute any longer and the finite squared amplitude is equal in 
all four regularization schemes. 

3 SCET approach to scheme dependence 

In Section 2 it has been shown that the regularization-scheme dependence of any massless 
QCD amplitude can be absorbed into a re-definition of the factor Z. Hence, it is important 
to study the scheme dependence of the anomalous dimension F governing the RG equation 
for the Z-factor. We work at NNLO, and at this order the anomalous dimension has a sum- 
over-dipoles structure. Thus, we need to compute the three relevant anomalous dimensions 
in Eq. (2.9), ycusp, 1 q and jg in the several schemes considered in this work, particularly in 
FDH (in DRED, also 7e is needed). In principle -jq and jg can be directly extracted from the 
IR divergences of the on-shell quark and gluon form factors computed in the three schemes. 
This approach [17, 18], which at first glance seems to be totally straightforward, turned 
out to hide highly non-trivial technical complications related to the UV renormalization 
procedure in schemes like fdh and bred. 

Here we show that the same y’s can be also extracted by combining the anomalous 
dimensions of the quark and gluon jet functions together with the anomalous dimensions of 
the corresponding soft functions (for Drell-Yan or Higgs production) defined through SCET 
operators. The soft and the jet functions can be computed with a standard diagrammatic 
procedure, and they are free of the renormalization difficulties that appear in the form 
factor calculations. This is an easier and more direct way to perform such a calculation. 
We have carried out this calculation at NNLO. In addition, the computation has also 
been carried out using the more traditional method to have an independent check of the 
results presented in this work and to show that the scheme dependence of these anomalous 
dimensions is universal and does not depend on the particular process analyzed. 

3.1 Outline of the method 

In the following we present the procedure for the direct calculation of the relevant anoma¬ 
lous dimensions in the four schemes via a SCET approach. The anomalous dimensions are 
obtained not from QCD scattering amplitudes but from soft and jet functions defined in 
SCET. Schematically, we get 


soft function ^ 

1 cusp; 

< h} ’ 

(3.1a) 

jet function => 

/cusp : 

^Rs 

’ ^9,9} ’ 

(3.1b) 


where h} Soverns the single-logarithmic evolution of the soft function for the case 

with an initial quark and an anti-quark (Drell-Yan) or two initial gluons (Higgs production), 
respectively. 7 ^® is defined similarly via the jet function. In bred, one has to distinguish 
the jet functions for D-dimensional gluons g and e-scalars g and the corresponding 7 °.^®° 
and The present discussion applies to these two cases in an analogous way. 

Thus, the cusp anomalous dimension 7 )[(fsp and its scheme dependence can be easily 
extracted independently either from the soft or the jet functions. The situation is slightly 


II 


more involved for the quark and the gluon anomalous dimensions where we need to exploit 
some known relations between anomalous dimensions to determine 7 ^® and 7 ^®. In the case 
of Drell-Yan and Higgs production, these relations hold as a consequence of the factorization 
of the cross section in the threshold region [43]. In particular one finds 

7^^,ov,h>=27,1,,+27{^), (3.2) 

where 7 ^® ^ is one half the coefficient of the 5(1 — x) term in the Altarelli-Parisi splitting 
functions and controls the parton distribution functions (PDFs) evolution. A similar re¬ 
lation involving the jet anomalous dimension instead of the soft anomalous dimension is 
found for DIS [44] 




'A9.9} 


o^,R-S 


(3.3) 


By combining Eq. (3.2) with Eq. (3.3) to eliminate the universal PDE anomalous dimension 
one obtains [43] 


= 7f 


7wj 


{DY, H} 


U^g} 


(3.4) 


The validity of Eq. (3.3) is a consequence of the factorization theorem for deep-inelastic 
scattering in the threshold region. The factorization proof is explicitly derived in [44] only 
for the quark current. Nevertheless by replacing the photon with a Higgs boson and after 
integrating out the heavy top loop, the factorization theorem for a gluon current follows 
in total analogy to the quark case. Indeed it can be explicitly checked that this relation 
holds both for the quark and gluon cases up to two-loop order by directly substituting the 
known expressions for the anomalous dimensions in cdr. 

Before we turn to the evaluation of the various anomalous dimensions we introduce 
some notation. As explained in Section 2.3 the anomalous dimensions in fdh and dred 
are equal, except for the appearance of the additional Y = 7 ^^*^°; see Eqs. (2.20) and 
(2.21). Likewise, the anomalous dimensions in CDR and HV are equal. Thus, we will drop 
the label RS whenever possible and denote fdh/dred quantities with a bar, schematically 


7 = 7 


= 7 


- _ „,FDH „,DRBD 

7 = 7 =7 


(3.5) 


In principle all perturbative expansions are carried out in terms of the five couplings {a}, 
as indicated in Eq. (2.1). However, for the results presented in this paper it is not necessary 
to distinguish the various Therefore, a coefficient in the perturbative expansion of 

the quantity X will have at most three labels, Ymnfcj indicating the power of Oe and 
Q! 4 e, respectively. Very often, the quantities do not depend on i.e. the last of the three 
indices is zero. In this case we often drop this label altogether and write the perturbative 
expansion with two labels only by setting Xmn = X^nO-^ 

We mention two special cases. Eirst, the f3 functions are defined with a negative sign. 




a . 

47r 


35 RS 
^mn 


(3.6) 


®In the CDR and HV schemes, all quantities of course only depend on Os. However, our notation will be 
adapted for the cases of fdh and dred, unless noted otherwise. 
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so the one-loop renormalization factors of as and ag in the various schemes are given by 




(3.7a) 

(3.7b) 


where the explicit form of the coefficients of the /3 functions are listed in Appendix B. 
Second, we also introduce an abbreviation for the cusp anomalous dimension multiplied 
with a colour factor, 


pRS _ RS _ 

cusp /cusp / y 


Uvr/ Uvr/ 


(3.8) 


where the colour factor Cr is either Cp or Ca, depending on the quantity under consid¬ 
eration. For brevity we omit the superscript cusp in the expansion coefficients of 

pRS 

cusp* 

3.2 Computation and scheme dependence of the soft functions and 'jw 

In this subsection we describe the calculation of the two-loop soft functions for Drell- 
Yan and Higgs production in momentum space and the extraction of the soft anomalous 
dimensions 'JWoy 7 Wh the different regularization schemes considered in this work. 
In the partonic threshold region, where the emitted gluons in the final state are soft, the 
Drell-Yan and Higgs production hard-scattering kernels factorize into the product of soft 
functions and hard functions. The factorization proof can be found in [26, 43]. The soft 
functions describe the real emission of soft gluons and contain singular distributions of the 
gluon energy while the hard functions depend on the virtual corrections and are regular 
functions of their variables. The soft matrix elements IF|nY.H}(^) arise in the cross section 
after the decoupling transformation which separates the soft and collinear sectors in the 
leading power SCET Lagrangian. 

The building blocks for the soft functions are the soft Wilson lines 


Si{x) = Pexp I igs / dsrii ■ Ag{x + 
\ J — OO 


srum 


(3.9) 


where A^{x) is a soft gluon held in SCET and n* = {n, h} (n^ = (1,0,0,1), = 

(1,0,0,-1) are light-like reference vectors in the direction of the two incoming partons). 
The path-ordering acts on the colour generators T“ in the representation appropriate for 
the zth held. For the conjugate quark helds one hnds T“ = — (t“)^ which turns into 
anti-path-ordering. The soft matrix elements IF{i 3 Y,H}(a^) are dehned in terms of a soft 
operator 

Os(x) = [SfiSn] (x), (3.10) 

as an expectation value of products of soft Wilson lines forming a closed Wilson loop 

VE{DY,H}(a:) = ^tr(0|f (Ot(x))r(O,(0)) jO), (3.11) 
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where dn = Nc for Drell-Yan and — 1 for Higgs production, T and T are the 

time-ordering and anti-time-ordering operators, respectively. 

Since the collinear and soft sectors no longer interact, it is worth noting that 
in Eq. (3.11) still contains the information about the colour and the direction of the initial 
quarks/gluons, but it is insensitive to the spin of the external particles due to the eikonal 
approximation. The soft function is defined as the Fourier transform of the soft matrix 
element 1 E{oy.h}(^) fo Eq. (3.11): 



dx^ 

Att 



0 ). 


(3.12) 


The Drell-Yan and Higgs production soft functions are closely related to each other; up to 
NNNLO they differ by Casimir scaling replacements [45]. At NNLO the situation is even 
simpler and the following replacement holds [46]: 


5H(a;) = + 0{al). (3.13) 

Thus, we directly compute the soft function for Drell-Yan and obtain the Higgs soft function 
by using Eq. (3.13). In the dred scheme the soft function for external e-scalars is also 
needed. Since soft gluon interactions are insensitive to the spinorial structure of the external 
particles, it turns out that the soft function for external e-scalars is the same as the one 
for external gluons. Therefore we will not discuss it further. 

In momentum space it is more convenient to rewrite the soft function in Eq. (3.12) as 
a squared amplitude by inserting a complete set of states 

(3.14) 

where Xg refers to a final state made of unobserved soft gluons carrying energy Exs • Foi* 
simplicity in Eq. (3.14) we drop the subscripts {DY, H}. To perform this calculation, we 
need not only the usual QCD Feynman rules but also the momentum-space Feynman rules 
for gluons emitted from Wilson lines up to O(a^). We report them in Figure 1. 

The 0{ag) [47] Drell-Yan soft functions in the CDR scheme have been originally calcu¬ 
lated in position space directly from the definition in Eq. (3.11). An exclusive soft function 
for Drell-Yan at 0{al) has been computed in [48]. The state of the art 0(a;f) soft func¬ 
tions for Higgs and Drell-Yan production have been computed very recently in a series of 
papers [45, 49, 50]. We also mention that related soft functions for thrust distribution and 
N-jettiness have been computed at 0{a‘l) in [51, 52] and [53] respectively. 

In order to study the higher-order corrections of the soft functions in the regularization 
schemes different from CDR we define expansion coefficients of the perturbative series as 


CRS 

'^bare 


(cj) — 6 {uj) + as (cj) Sfff (lu) + a‘l{u}) S20 (w) -|- ... , 


(3.15) 


where we have introduced the superscript rs to indicate the scheme dependence. In the 
above equation we have introduced 


as{uj) = e 



^bare 


47r 


J (47r) 


(3.16) 
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Figure 1. Feynman rules for the emission of one and two gluons from a Wilson line. Figure taken 
from [26]. 


and expressed the bare coupling in terms of the renormalized coupling as = Q!s(Ai) 

in the MS scheme. Note that as{uj) and are actually scheme independent, but if 

expressed in terms of the MS coupling as{n) depend on the scheme-dependent renormal¬ 
ization factor The all-order bare soft function in Eq. (3.15) is independent of the 

renormalization scale Up to NNLO the soft function depends only on Og and not on Oe 
or Q;4e. 

At NLO only two diagrams contribute to the soft functions; they describe the real 
emission of one soft gluon from the Wilson lines. At NLO the bare soft function turns out 
to be scheme independent, 


SioH 


8 e®T'^r(-6) 

u^T{l-2e) ■ 


(3.17) 


As a result, the soft anomalous dimensions must be scheme independent, too. This re¬ 
produces the well-known fact that 7cusp is scheme independent at NLO, and it implies 
7 ^^® = 0 in all RS. The reason is that for the fdh and dred schemes there are no ad¬ 
ditional diagrams involving e-scalars compared to CDR and HV. This is a consequence of 
the fact that dot products of a e-scalar field A with the vectors n, n are vanishing, i.e. 
n ■ A = n ■ A = 0. It follows that soft e-scalars cannot be emitted from the Wilson lines. 
This explains in a direct way the result [4] that the scheme dependence of general NLO 
amplitudes is contained in the parton anomalous dimensions. 

At NNLO the situation is more involved; diagrams with two real soft emissions and 
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virtual diagrams with one real soft emission are present. The soft functions and soft 
anomalous dimensions at NNLO have a scheme dependence, which originates from the e- 
scalar cut bubble contributing to the second diagram in Figure 2. The grey blob represents 
the quark, gluon, ghosts and e-scalar contributions. The latter is present only in fdh 
and DRED. After calculating the non-vanishing integrals using the techniques described in 
[54, 55] and summing all the contributions we obtain the NNLO coefficient in Eq. (3.15) 
in fdh/dred. 


52o(w) = -Cr [Ca Sa + NfTr Sf + Cr Sr] 


OJ 


with 


^ -ur + 


44 2N, 




1 / l6iVe 47r^ _ 268 \ 


e V 9 


9 / 


Tvr^iV, 104A', 
+ 


+ 56C3 + 


1547r^ 

9 


1616 


+ - 


1241V,C3 567r2iVe 


640A^e 2728C3 47r^ 

H-^ H-^-h 

81 9 9 


9 27 

9387r2 9712 \ 

^7 sF ) 


e + 0 (e 2 ) 


- 16 8 0 567r2 448 


992(3 2624 2807r2 


9 


81 


27 


32 1127r2 1984(3 4Fe 


5k - - ^ ^ + O (e^), 


27 


e + 0 (e 2 ) 


(3.18) 


(3.19a) 

(3.19b) 

(3.19c) 


where Cr = Cr for Drell-Yan and Cr = Ca for Higgs production. 

We now turn to the determination of the soft and cusp anomalous dimension from 
the soft function. In order to do this we need to discuss the singularities of the soft 
function that remain after coupling renormalization. From the point of view of ordinary 
QCD computations, these remaining singularities are closely related to IR singularities. 
However, from the SCET point of view they simply correspond to UV singularities and are 
to be removed by renormalization within the effective theory. For convenience this is done 
in Laplace space by introducing the Laplace transformed soft function as 


F®(/^) = ^“du;exp (-^) 5«®(a;), 
where the integral transform can be easily carried out by using the relation 



doj exp {—boj) oj 


1—ne 


r(-ne) 6 ”^. 


(3.20) 


(3.21) 


The remaining UV divergences of the soft function can be subtracted multiplicatively. 




(3.22) 


Like in the case of general amplitudes in Eq. (2.6) and Eq. (2.7), the RGE 


d 

d\n^ 








(3.23) 
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Figure 2. Selected non-zero Feynman diagrams contributing to the one-loop and two-loop soft 
functions. A complete list of diagrams can be found in [47]. Double lines indicate the direction of 
Wilson lines while the red vertical cut indicates on-shell partons. The scheme dependence originates 
from the diagram D 2 ■ Diagrams D 2 , D 3 and D 4 represent double real soft emissions while diagram 
D 5 represents a single virtual-real emission. 


holds, and the corresponding anomalous dimension has a structure similar to Eq. (2.9), 

^^) = [-4r[(4p Lk - 27 ]^] > (3-24) 

which is derived from the RG invariance of the cross sections in the threshold region in 
analogy to the CDR case in Ref. [43]. In Eq. (3.24) we have defined = ln(K//i) and 

Cr = Cp for Drell-Yan and Cr = Ca for Higgs production. Comparison of the previous 
two equations yields an expression for the fdh renormalization factor Zs{k, fi) = 
in terms of the soft and cusp anomalous dimensions. This expression has the same structure 
as Eq. (2.12), but can be written in a simpler form because up to NNLO the soft function 
does not depend on Oe and Q: 4 e: 


In Z. = 


= 

Vdvr 


+ 

+ 0{a 


dvr/ 
3\ 


-^ + t{2r.„L. + 4f!;) 

3 / 320^10 


(3.25) 


4g3 


2e2 


{2TiQLf^ + + — {2V2 oLk + 7 ^) 


By requiring that the renormalization factor Zg in Eq. (3.25) minimally subtracts all of the 
divergences of the bare soft function (in fdh, treating as an independent multiplicity). 
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we extract the expressions for the anomalous dimensions in the fdh scheme 


r — 

-*■ cusp — 


+ 


dvr 

Vdvr 


Cr(4) 


Cr 


Ca 


268 


+ 0{a^). 


3^ 
808 11 


- f THiVf - N.fc^ 

52 


TT' 


+ O(a^) 
2 , 


Ca - ^ —vr^ + 28C3 + 


27 


9 


27 
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(3.26a) 

- 1a)1 

27 9 / 

(3.26b) 


The fact that f cusp = C/j7cusp; with the known expression of the cusp anomalous dimension 
in the FDH scheme, 7cusp) is a consistency check of the method, is a new result. The 
corresponding expressions in cdr/hv can be obtained by simply using the appropriate (3 
functions and anomalous dimensions in Eq. (3.25) and by setting = 0 in Eq. (3.26). 
They are consistent with the literature [43]. 

Finally we remark that in analogy to Eq. (2.5) we can obtain a finite and scheme 
independent soft function Sfin through 


Sfin(«^,7t) = , lim (3.27) 

(iV)E^O 

The explicit expression for sgn is given in Eq. (A.2) of Appendix A. 


3.3 Computation and scheme dependence of the quark jet function and yjg 

The quark jet function has been calculated at NNLO in cdr [28]. Referring to [28] for 
more details, we describe here the corresponding calculation in fdh (which is identical to 
the one in dred, but for simplicity we will only refer to fdh in the present subsection). 
The jet function is given in terms of the hard-collinear quark propagator 

= j (fxe"P'^{()\T{xhc{x)Xhc{^)}\^) 

= j , (3.28) 

with Wilson lines ^ 

W{x) = V exp (aqs / dsn ■ A{x + sn) \ , (3.29) 

^ J —OO ' 

where = Aa t°‘. The field Xhc{x) is the gauge-invariant (under both soft and hard- 
collinear gauge transformations) effective-theory field for a massless quark after a decou¬ 
pling transformation has been applied, which removes the interactions of soft gluons with 
hard-collinear fields in the leading-power SCET Lagrangian. As shown in Eq. (3.28), we 
can rewrite the propagator in terms of standard QCD fields. 

The hard-collinear quark propagator as defined in Eq. (3.28) is scheme dependent. 
The fields Xhc and -0 on the r.h.s. of Eq. (3.28) are Heisenberg fields, so applying the 
usual perturbative expansion results in loop diagrams contributing to the propagator. The 
scheme dependence is related to UV singularities of such diagrams. Examples of two-loop 
diagrams are shown in Figure 3. In fdh the computation is similar to the CDR scheme. 
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Figure 3. Examples of two-loop diagrams contributing to the quark jet function. Gluons emitted 
from the crossed circles originate from the Wilson lines. Diagram (a) contributes in CDR and fdh, 
whereas diagram (b) with two e-scalars contributes only in fdh. 


However there are additional diagrams, which include the e-scalars and also depend on 
the coupling a^- An example of a two-loop diagram needed for the jet function in fdh 
(and not present in the CDR scheme) is shown in Figure 3 {h). Since h is a ZD-dimensional 
vector, there are no e-scalars originating from the Wilson lines. Indeed, the scalar product 
in Eq. (3.29) will vanish in the case of the e-scalar. 

The jet function J^®(p^) is the discontinuity of the propagator, i.e. 


^ vr 


ijr 


(p 2 


(3.30) 


To highlight the similarities with the discussion in Section 2 and the soft function it is 
convenient to work in Laplace space, so we define jf^(Q^), the Laplace transform of the jet 
function as 

(Q2) = ^ dp^ exp 

The analogous equation in the case of the soft function is Eq. (3.20). 

To compute the propagator in the fdh scheme, Sqip^), the diagrams have been gen¬ 
erated with QGRAF [56] and the colour algebra has been done with ColorMath [57]. For 
the reduction of the integrals Reduze 2 [58] has been used. The master integrals needed 
for the FDH jet function are the same as for the CDR scheme. After taking the imaginary 
part and performing the Laplace transform, the bare quark jet function at NNLO in fdh 
is obtained as 


Jgbare(Q^) — 1 + Os(Q^) Cp g -3-^ e(l4 - - -8^3)^ 

+ ae(Q^) Cp (^ — — — 1 + — 2 + 

+ a2(Q2) + CpCAllt + CpTpNpjg) 

+ aliQ^) + CpCAlSt + CpTpNpj^i) 

+ as(Q^)oe(Q^) Jii ~\~CpCaJi{ ^-b O(a^). 


In analogy to Eq. (3.16) we have defined 


as{Q‘^) = e-^^^{ApY 



^bare 


47r 


\ ^OLs^S 

QV (47r) 


(3.32) 


(3.33) 
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with an analogous equation for Og. The explicit expression for the two-loop coefficients are 
given in Appendix A. Note that jqbare(Q^) is independent of 

The renormalization procedure in any regularization scheme can easily be general¬ 
ized from the corresponding procedure in CDR [28]. A renormalization factor 
absorbing the UV divergences of the bare jet function is introduced such that 




(3.34) 


is finite. This equation is analogous to Eqs. (2.4) and (3.22). Requiring minimal subtraction 
with as an independent multiplicity determines the explicit form of ZJ^ {Q‘^, fj.) uniquely 
in terms of the bare quark jet function Jgbare(Q^)- In principle, Zj® depends on all couplings 
{a}. However, in fdh, up to NNLO there is no dependence on 04 ^. 

To relate to the cusp anomalous dimension and the quark jet anoma¬ 

lous dimension 7 j® we follow the same procedure as for the soft anomalous dimension. We 
compare the RGE of the quark jet function in the form 


d 




dZf{Q\^i) 


dln/i dln/r 

to the RGE written in terms of r[(®gp and 7 j®, 


zXiQ\d)) 


-1 




(3.35) 


d 

dln/r 





RS 

cusp 




(3.36) 


This relation is analogous to Eqs. (2.6) and (3.24); we have used Lq = ln(Q^///^) and 
r^isp = Cplmsp- With the help of Eqs. (3.35) and (3.36) we can express Zj^ in terms 
of the FDH anomalous dimensions. Up to NNLO, the expression for In Zj^ has the same 
structure as Eqs. (2.12) and (3.25). We write it explicitly, using that up to NNLO only the 
two couplings as and ceg appear: 


1 7 


1^10 I I /"P r I —jq 

- 3 - + - (Tio Lq + 7 ^q 


-k 


a^ 
dvr . 


TOI ,1 /f=, j . -Jq 
- 3 - + - ( Toi Lq + 7o^ 


[ 3 (/3|orio-k/^loToi) /3, 


4^3 
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4g3 
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Toi Lq + 7qi 

(3.37) 


On the one hand this formula gives strong consistency checks. It allows for an independent 
extraction of the cusp anomalous dimension and the coefficients of the /3 functions of ag 
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and «£ in the fdh scheme. These coefficients agree with the well-known results in the 
literature [17, 18]. 

On the other hand, comparing Eq. (3.37), in particular the 1/e pole, to the explicit 
result for the bare quark jet function allows to read off the anomalous dimension 7 j^. We 
obtain the following explicit expression in the fdh scheme: 




7^= ^ (-3Cf)+ P -^Cf 
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dvr/ 

dvr 


Ue 
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(3.38) 


Using this expression together with Eqs. (3.d) and (3.26b) the quark anomalous dimension 
in the fdh scheme, 7 ^ can be found. Thus the computation of the soft and quark jet 
functions provides an alternative determination of 7 ^. The result agrees with previous 
determinations [17, 18] and is listed in Appendix B for completeness. Of course, setting 
W = 0 only the pure ag terms survive and the well known results in the cdr/hv scheme 
are recovered. 

This is also true for the quark jet function as a whole. In analogy to Eq. (3.27) we can 
define 


Jgfin(Q^/^) = ^ lim (3-39) 

{N)e^0 

SO the finite quark jet function is scheme independent and can be obtained using any of 
the regularization schemes. The explicit result is given in Appendix A. 


3.4 Computation and scheme dependence of the gluon jet function and jjg 


The discussion of the previous subsection can be readily adapted to the gluon case. We 
closely follow Ref. [29], where the gluon jet function Jg{p^) has been calculated at NNLO 
in CDR. The starting point is the gauge-invariant field related to the collinear gluon 
field Ac{x) through 

A^{x) = A^^{x)ta = W\x)[iDi^W{x)]. (3.40) 


The treatment of this vector field depends on the regularization scheme; we will give the 
details below. In all schemes the field A^ satisfies n • ^ = 0; hence it can be decomposed 
as A^ = A^^ -|- (n • JC)fAI 2 and the leading term is A'^. The gluon jet propagator Jg{p^) 
is then defined as 

i-gT) = / d^xe^P^{0\T{AT{x)A’Tm\0). (3.41) 

For the calculation of it is actually more convenient to use an equivalent definition 

in terms of the time-ordered product of the full fields A^, 




— gfj-u + 




n ■ p 


(n • pY^ 


(3.42) 
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Figure 4. Sample two-loop diagrams contributing to the gluon jet function. Diagram (a) is 


present both in CDR and fdh, diagram (b) including an e-scalar contributes only in fdh. 


/ 


= / d^xe^P^0\T{A^^{x)Alm\0) 


and then extract using a projection. The gluon jet function is the disconti¬ 
nuity of the leading part of the propagator, more precisely = Im[i The 


function is related to power-suppressed terms and will not be considered any further 
in this paper. 

As in the case of the quark jet function, after decoupling of the soft fields, the collinear 
Lagrangian is equivalent to the QCD Lagrangian. Exploiting the gauge invariance of 
we work in the light-cone gauge h • A = 0. This is particularly convenient as in this gauge 
W(x) = 1 and, therefore, no diagrams with additional emission of gluons from the Wilson 
lines have to be considered. Therefore, for the calculation of only standard QCD 
Feynman rules are required. Of course, ghost loops are also absent in this gauge. 

Now we give details on the regularization scheme dependence. Typical examples of 
two-loop diagrams contributing to are shown in Figure 4. In cdr all gluons are D- 
dimensional gluons g and no e-scalar diagrams are present. Correspondingly, the metric 
tensor in Eq. (3.42) is in CDR. In HV and fdh the external gluons are understood to be 
strictly 4-dimensional. Thus, the gluons attached to the Wilson lines in Figure 4 are to be 
interpreted as g, and the metric tensor in Eq. (3.42) is g^'^ in these schemes. Furthermore, 
in FDH internal gluons are treated as g and hence are decomposed into g and g, as indicated 
in the left and right panel of Figure 4. In dred the definitions of the present subsection 
apply to external D-dimensional gluons g. For these, the calculation and the result are the 
same as the corresponding fdh calculation, see Eq. (2.20). Hence for simplicity we will 
only refer to FDH in the remainder of the subsection. 

After an explicit calculation of the diagrams in fdh, taking the imaginary part and 
performing the Laplace transform, we obtain for the bare gluon jet function in fdh 



+ a? (^^1 + CaNfTrJI^^^ + + iV|r| 

+ o,ae + CFNFTjijfi^^^ + 0{a^). 


(3.43) 


22 












The explicit results of the two-loop coefficients are given in Appendix A. In the limit ^ 0 
all terms proportional to ag vanish and we obtain the results in CDR, in agreement with 
Ref. [29]. 

The renormalization procedure is the same as for the quark jet function. In Laplace 
space, the renormalized gluon jet function in the fdh scheme is obtained by multiplying 
Eq. (3.43) by a factor Zj^. This factor is the same as in Eq. (3.37) apart from the re¬ 
placement —>■ 7 ]^® and . After renormalization of the coupling, all 

divergences of the bare gluon jet function have to be absorbed by This allows 

to determine the anomalous dimension of the gluon jet in the fdh scheme as 
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(3.44) 


Of course, it is again also possible to extract the cusp anomalous dimension as well as the 
j3 functions of as and ag from Zjg{Q‘^,iJ,). The fact that we obtain again the same results 
for these quantities is a strong consistency check on the procedure. 

From 7 jg we can determine 7 ^ with the help of Eq. (3.4). The result is in agreement 
with previous determinations [17, 18] and is listed in Appendix B for completeness, but 
the present procedure provides a more direct alternative determination of 7 ^. 

Finally, as for the soft and quark jet function, we can obtain a finite and scheme 
independent gluon jet function as 


Jafin(<3^1i) = ^ lim JgLb(<3^/^)• (3-45) 

(Af)e^O 

For completeness the explicit result is listed in Appendix A. 


3.5 Computation of the e-scalar jet function, jjg and result for 7 c in DRED 


In DRED processes with external e-scalars need to be considered. The discussion of Sec¬ 
tion 3.1 applies analogously, and we can determine the anomalous dimension of e-scalars 
from an equation like Eq. (3.4), 


DRED 

ig 


= = 7J. - 


— ,DRED 
'We 


(3.46) 


As mentioned in Section 3.2 the soft function is the same as for external gluons, hence 
7 t^^° = 7 vi/) from Eq. (3.26b). For 7 j^ an e-scalar jet function is needed. Such an object 
can be dehned and computed in close analogy to the calculation of the gluon jet function, 
with the difference that now the time-ordered product of two fields = g^j.uA'^ has to be 
considered. In light-cone gauge these fields reduce to the e-scalar field A^. Starting from 
the propagator Se{p^) = given by 


i-~9g.) Up") = I d^xAP^0\T{A<;:,ix)Alm\0) (3.47) 
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Figure 5. Sample two-loop diagrams contributing to the e-scalar jet function both. Diagram (a) 
is proportional to OsOfe whereas diagram (b) is ^ 


the e-scalar jet function is obtained as J^ip^) = \va{i 

Two examples of diagrams contributing (in light-cone gauge) at two-loop order are 
shown in Figure 5. A new feature is the appearance of the quartic coupling We do 
not need to distinguish the three different a 4 e since the quartic coupling only appears at 
the two-loop level and hence the associated renormalization constants and /3 functions do 
not appear. The only non-vanishing diagram ~ a\^ is depicted in Figure 5 b. 

Performing a computation analogous to previous cases, the bare two-loop e-scalar jet 
function in Laplace space is found to be 

Jebare(Q^) = 1 + 05 Ca + - + 8 -^-h e(l6 - - - 

+ Oe — — —4-|-€( — 8 -|- 

+ «5 (^^1 + CANpTnJ^i,^^) 

+ al NfTr [Ca + Cf Jo20^ + ^fTr 


I 2 /o 2 —-A- 

+ ®4e J002 


+ a,Oe NfTr (Ca Jl[o^ + Cf + 0{a^). (3.48) 

Due to the presence of a 4 e, the various coefficients have now three labels, with the last one 
indicating the power of a 4 e. The explicit NNLO expressions are given in Appendix A. 

Once more, the UV divergences of the bare jet function are absorbed by a renormal¬ 
ization factor {Q'^, fi), which has a structure similar to Eq. (2.12) or Eqs. (3.25) and 

(3.37). In fact, it can be written as Eq. (2.18), 


In^DRED ^ 



(3.49) 
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with the identification 


p.DRED ^ -4^CAicus^, 


= 2 Ca Tcusp Lq + 2 7j, . (3.50) 
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We refrain from using the explicit form of Eq. (3.37) since the dependence on 04 ^ leads 
to a proliferation of similar terms. The only simplification used is the identification of the 
couplings which is possible since the explicit results show that these couplings appear 
not at one-loop but only in the genuine two-loop coefficients. 

By comparing with the explicit result for the e-scalar jet function we determine the 
renormalization factor using minimal subtraction and extract from this the anomalous 
dimension of the e-scalar jet as 
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(3.51) 


Combining this result as prescribed by Eq. (3.d6) with the soft anomalous dimension, which 
has only contributions, we find the e-scalar anomalous dimension 
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(3.52) 


As discussed in Section 2.3, 7 ^ is needed to relate two-loop matrix elements computed in 
DRED to those computed in other schemes such as fdh. With this new result all anomalous 
dimensions are known at the two-loop level in all four schemes. 


4 Alternative determination of 7 ^ from the e-scalar form factor 

Apart from the new approach of extracting the IR anomalous dimension of the e-scalar, % 
defined in Eq. (2.21) from the e-scalar jet and soft functions, it is also possible to obtain 
this quantity in the more traditional way, by comparing the generic infrared factorization 
formula with a specihc amplitude for a process containing external e-scalars. This pro¬ 
cedure is analogous to the determination of 7 ^ and 7 ^ in Ref. [18]. We now describe the 
determination of % via a process with two external e-scalars, the e-scalar form factor, which 
has been calculated recently in Ref. [31] up to the two-loop level. 
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According to Eq. (2.18) the one-loop infrared divergences in the dred scheme are 
described by 


InZ^L 
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2e2 ^ e 



^010 7^ 

2e2 ^ e 





^001 7^ 

2e2 ^ e 


(4.1) 


Here the relations T'ijk = —2Tijk = —2CA^lJk^ ^ijk = ‘^itjk have been used. The 
notation with three indices for a common a 4 e coupling and for dropping the superscript 
“cusp” has been explained in Section 3.1. Eq. (4.1) can now be compared with the corre¬ 
sponding IR divergent one-loop result of the UV renormalized e-scalar form factor given in 
Ref. [31], where Tr = ^ and = —si 2 has been used: 
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(4.2) 


The ^-pole of this one-loop form factor confirms the previous hnding that the one-loop cusp 
anomalous dimension is a process-independent quantity that has only one non-vanishing 
component Tioo = 4 Ca = 7ioo^ Ca- On the other hand, the 4-poles in Eq. (4.2) are directly 
correlated with the components of the anomalous dimension The values obtained here 
agree with the results from the previous section. 

The appropriate two-loop prediction for In Z2^ could be given in a completely general 
form, as in Eqs. (2.18) and (3.49), in which it would allow to read off once again even the 
one-loop /3 functions. Here, however, we give the prediction in a more specific form, where 
we already use the knowledge that several one-loop coefficients are zero. Considering only 
non-vanishing components of one-loop anomalous dimensions and (3 functions yields for 
the infrared divergence structure at the two-loop level: 
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(4.3) 


Thanks to the simple colour and momentum structure of the form factor, this has to 
correspond directly to the divergence structure of the combination .F2L _ i||^^iL^ 2 ^ 
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Ref. [18]. Inserting the results for the form factor of Ref. [31] yields 
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(4.4) 


Again, the ^-poles allow to read off the components of the anomalous dimension of the 
e-scalar %. The values found here agree with the results from the previous section, see 
Eq. (3.52). Since the remaining divergence structure is governed by one-loop anomalous 
dimensions, the process-independent components of the cusp anomalous dimension and 
previously known (5 coefficients, this is further evidence for the validity of the results 
obtained in Section 3.5. With this result, and the results of the previous sections and 
Ref. [18], all two-loop anomalous dimensions 7 * in all RS have been determined both in the 
SCET approach and from form factors. 


5 Cross check with explicit processes 

The results of the previous sections allow us to predict the differences between UV renor¬ 
malized virtual two-loop amplitudes squared, as defined in Eq. (2.3), computed in different 
regularization schemes. In this section we will make these transition rules more explicit 
and will check them with explicit examples. 

The following discussions will also shed more light on the role of the various couplings 
cts, Oe and a 4 e,j. In the practical computation of the genuine two-loop diagrams it is no 
problem to set these couplings equal from the beginning. In the process of UV renormal¬ 
ization, i.e. in lower-order diagrams with counterterm insertions, the bare couplings and 
the associated renormalization constants appear. It is unavoidable to keep these distinct, 
regardless whether fdh or dred is used. Once renormalization has been performed, it is 
possible to set the renormalized couplings equal and to identify W and 2e. Likewise, the 
derivation of the IR subtraction formulas and the transition rules requires the couplings to 
be treated independently, but in the end the transition rules can be easily written down 
for the special case of equal couplings. 

We will consider the transition rules fdh ■<-)■ HV, as well as fdh •<->■ dred. To make con¬ 
nection to the scheme that is used most often, cdr, we remind the reader of the discussion 
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in Section 2.2. The only difference in the squared matrix element between HV and CDR is 
due to the use of different metric tensors for the polarization sum of external gluons. All 
anomalous dimensions are the same in the two schemes. 

5.1 Transition between FDH and HV 

Since external gluons are treated in the same way in fdh and HV, we can actually relate 
directly virtual amplitudes and do not need to work with squared amplitudes. The finite 
remainders of the scattering amplitudes are scheme independent. More precisely 

= limZ"^(e,{p},/r)|M(e, {p})) 

= lim Z~^{e,N^,{p},fi)\A{e,N^,{p})), (5.1) 

(Af)e-i.O 

where |M) = and Z = denote quantities in the HV scheme and |.A) = and 

Z = are the corresponding quantities in the fdh scheme. Suppressing the arguments 
of the amplitudes, setting = 2e and writing Z“^ = 1 + dZ in both schemes, we can 
rewrite this equation as 


\A) + 5Z\A) = \A) + 5l\A) + 0{e). (5.2) 

If the expansion coefficients dZ are known to 0{a"') and the amplitudes |M) are known to 
this equation allows to obtain a relation between the C>(a”) amplitudes computed 
in HV and fdh, up to 0{e) terms. We now give the explicit results up to the two-loop level. 

The tree-level amplitudes in the two schemes are the same |Mo) = |.Ao). At one-loop 
we can relate the 0(as) and 0{ae) corrections in the fdh scheme, denoted by |.Aio) and 
IMoi) respectively, to |Mi), the 0(as) corrections in the HV scheme 

|Moi) = —5Zqi\Aq) + 0(e ), (5.3a) 

lAo) - l-Ai) = (5Zi - <iZio)|A) + 0(e}. (5.3b) 

In the above equation we have also introduced the expansion coefficients 5Zm and dZ^n of 
Z~^ = 1 6Z in the HV and fdh scheme, respectively. Substituting in the last equations 
the explicit expressions of these expansion coefficients, the explicit form of the differences 
for a process with external massless quarks and external gluons read 

lAi) = ^1^1 A) + 0{e) = + 0{e ), (5.4a) 

l^o) - I A) = + 0{e) = #5^|Mo) + 0{e ), (5.4b) 

which agrees with the results in [4, 16]. In Eq. (5.4) and what follows we use the notation 
(see footnote in Section 3.1) 7^0 = for the anomalous dimensions (and the /3-functions) 
in the hv scheme. Since in the hv scheme the anomalous dimensions depend only on ag but 
not on «£ the second label is always zero. Of course, this is not the case in the corresponding 
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quantities in the fdh scheme, 7 mn- To obtain Eq. (5.4) we have used and 

cusp -cusp 

7io = 7io • 

Moving to the two-loop level the corresponding equations are 


I-T 02 ) — —<5Zoi|Moi) — 5 Zo2|^o) + C?(e), (5.5a) 

l^o) - I-T 2 ) = 5Zi\Ai) - 5Zio|Ao) + (- 5 Z 2 - 5Z2o)|A) + 0{e) , (5.5b) 

|Mii) = - JZiolAi) - <5Zn|Mo) + 0{e) . (5.5c) 

The expressions given in (5.5a), (5.5b) and (5.5c) allow one to move from fdh to HV (and 
vice versa) for any process with external gluons and external massless quarks in 
QCD up to two-loop order. Exploiting and = 7io'^’’ we obtain 
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where we have defined 


\Aif) = \A,o)-\A,), 


(5.7a) 


1-4' 


finx 


= lim 

e-i-0 


5'Li\Ao) + l^i) 


= lim 


(5Zio|^o) + |-4io) 


(5.7b) 


|^f“) is the NLO approximation to |^fin) and, thus, a finite and scheme independent 
quantity. The one-loop quantities |^io®) and |^oi) have to be known up to O(e^) terms. 

We remark that Eq. (5.5a) allows to obtain the 0(ag) contribution of a two-loop 
amplitude in fdh up to 0(e} terms directly from the tree-level amplitude. This is due to 
the fact that ~ W ~ e and hence the coefficient multiplying |yloi) in Eq- (5.6a) is 
finite. Therefore, we can use Eq. (5.3a) and with the explicit expressions of the anomalous 
dimensions we get 


\Ao2) = Cf 
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6 NfTf) 


l-^o) + C>(e). 

(5.8) 


Eor a process with no external quarks, = 0 there are no 0{a‘^) terms at NNLO, as can 
easily be confirmed on a diagrammatic level. 

As mentioned several times, once the UV renormalization has been carried out, there 
is no need any longer to distinguish between the different couplings. After setting = OLg 
the full difference is given by 
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Figure 6. Examples of two-loop (left panel) and one-loop counterterm (right panel) diagrams for 
99 99 (top panel) and qq —> gg (bottom panel). Black vertices denote couplings gs whereas 

white vertices denote couplings ge, and crosses denote counterterm insertions. For gg —>■ gg at 
one-loop, there are no contributions with couplings ge- The order is given relative to the Born term 

lA) - 0{a,). 


where we have introduced the notation 

1-^2) = \A 2 o) + |.4o 2) + |-4ii), (5.10a) 

= l“4io) + |-4oi) — l^i). (5.10b) 

5.2 NNLO 2 —> 2 amplitudes in HV and FDH in massless QCD 

As an example for the transition rules derived in the previous subsection, we consider the 
two-loop amplitudes gg —>■ gg and qq —)■ gg for massless quarks. Initially the interference 
of these two-loop amplitudes with the tree-level amplitudes was calculated in cdr [59, 60]. 
Later the helicity amplitudes were computed and explicit results in the HV and fdh scheme 
were given [61, 62]. However, for the computation and the UV renormalization procedure in 
the FDH scheme, no distinction between a<j and Oe (and was made. For the process 

99 99 this is of no consequence, but for qq gg this will lead to an incorrect UV 

renormalization. As shown in Refs. [6, 7, 17] this leads to incorrect finite terms which 
violate unitarity. For our purposes it also matters because an incorrectly renormalized 
amplitude cannot be consistent with the IR structure and transition rules discussed above. 

Hence, in order to check the validity of the transition rules we first need to correct 
the renormalization of the qq —>■ gg result of Ref. [62]. Figure 6 shows diagrams which 
illustrate the problem. The left panels show genuine two-loop diagrams to gg —>■ gg and 
qq —)■ gg. One of them depends on cie, but setting = «« in these two-loop diagrams 
causes no problem. However, the diagrams have subdivergences, which should be cancelled 
by suitable counterterm diagrams, such as the ones in the right panels. The first of these 
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counterterm diagrams depends on the one-loop renormalization constant but the 

second one depends on which differs by a divergent amount. If, as in Ref. [62], this 

renormalization constant is effectively replaced by dZ^^, the subdivergence is not properly 
subtracted, and the final result will not be correct. 

The correct renormalization procedure requires to compute the lower-order amplitudes 
for individual couplings. At tree-level, the amplitudes |.Ao) for both processes are propor¬ 
tional to as and hence are correctly renormalized by multiplying with . At the one-loop 
level, the amplitudes receive contributions of 0{as) or 0{ae) relative to tree-level. The 
latter contribution |^oi) must be renormalized by multiplication with Za^Za^. 

The difference between the two processes gg gg and qq gg is that for the former 
process, j^oi) happens to vanish. This is the reason why for this process the identification 
as = Oe causes no problem. In order to restore the correct renormalization for the latter 
process, we have computed the 0{as Oe) contribution to the one-loop amplitudes. We have 
then renormalized this contribution using Z^^Z^^ and add the resulting NNLO term to the 
explicit results of Ref. [62] . We also subtracted the corresponding terms obtained with the 
renormalization factor Z\^ that had been applied in Ref. [62]. 

We have compared the difference between the fdh and HV amplitudes for both pro¬ 
cesses with the prediction given by Eq. (5.9) and have found full agreement. This is a 
further non-trivial confirmation that our treatment of the scheme dependence is process 
independent and applicable at least to NNLO. It is also an independent verification of the 
correctness of the anomalous dimensions in fdh. 


5.3 Transition between FDH and DRED 


The transition rules between dred and fdh can be derived similarly but are more involved. 
To illustrate their structure let us first consider a process with a single external gluon. The 
explicit calculation of the UV renormalized matrix element in dred yields that 

can be written as 

M°^^°{g)=M°^^°{g)+M°^^°{g) = 2Re{Al\A^) + 2Re{A^o\-^^), (5.11) 

where we have introduced the shorthand notation A^ = .4.°^™(5) and A^ = .4.°^™(5) etc, 
and suppressed other arguments compared to Section 2.3. We would like to find a relation 
between and the corresponding result in fdh, 

M™«(<7) = 2Re {Al\A^) = 2Re (^S°"(5)l-4™"(5)) • (5.12) 

To do so, we start from the equality of the IR subtracted amplitudes computed in dred 
and FDH, written with a similar shorthand notation for the Z-factors as 

KI{Z>)-‘|» + (.45|(Z-)-‘|.4-> = + 0(e), (5.13) 

where we have set W = 2e. Writing = 1 -|- dZ, where dZ denote the perturbatively 
expanded higher-order terms we obtain an equation analogous to (5.2), 


M°^^^{g)+2R(i{Al\6Z3\A9) + 2Re{Al\6Z^\A^) 
= M^^^{g) + 2Re{Al\5Z^\A9) + 0{e). 


(5.14) 
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If the expansion coefficients dZ are known to 0{a^) and the amplitudes |^) are known 
to Eq. (5.14) allows to obtain a relation between the 0{a^) squared matrix 

element computed in dred and fdh, up to 0(e) terms. For this relation, the knowledge 
of TJ" = is required, even though Eq. (5.13) is still correct if the second term on 

the l.h.s. containing Tf is dropped. 

As a concrete example we consider the process H g g m. fdh and dred and work 
out the transition rules between the two schemes for the UV renormalized two-loop squared 
amplitudes. For simplicity we also set tte = ot^e = «s- 

As we have i^g = 2 external gluons, in dred the squared matrix element is to be 
written as a sum over 2^^ = 4 terms. However, in this particular case two of these terms 
vanish to all orders, resulting in 

M°^'^^(g,g)=M(g,g) + M(g,g). (5.15) 

Writing explicitly the equality of the subtracted matrix elements in fdh and dred we get 



+ 0{e) + 0(a^). (5.16) 

In Eq. (5.16) we have introduced a compact notation for the perturbative coefficients of 
the amplitudes and Z“^ in dred: = {• A . 2 ( g , g )) and 

Z-i (g, g) = l + 5Zf (^) + 6Zf + 0(a3), (5.17) 

with analogous expressions for other partonic processes. Comparing the order Og terms 
yields 

Mr^^(g,g) - Mr^(g,g) = Mr^^(g,g) ~ 

= + 0(e). (5.18) 

This one-loop transition rule is in agreement^ with Ref. [4]. To make this agreement more 
explicit we write the transition in a more general way as 

(5.19) 

Note that the difference is finite, since the tree-level matrix element squared on the r.h.s. 
of Eq. (5.18) or Eq. (5.19) are of 0(e). 

“^Note that in Ref. [4] a different convention for the y’s has been used. 
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In order to write the scheme difference at NNLO we introduce a similar short-hand 
notation for the squared matrix elements as for the amplitudes, denoting the full tree-level 
and one-loop contribution for the H ^ gg process by A4q'^ = Moiff, 9) and = M.i{g, g), 
respectively. The difference can then be written as 


M^^^°ig,g) - Mf^^{g,g) = 7 ^ (Toio + 7mo “ 7ioo) 


2e2 I 


2e3 

■^ 0 "(^0207010 + ^iio7oio + ^ 200(7100 - 7ioo) + (7oio + 7ioo)^ “ (7ioo)^) 
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2 n 


+ CAMriZ^ - C^A^“7ioo"7mo 
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2 A^f (7010 + 7ioo “ 7ioo) + -^ 0^(7002 + 7o 20 + 7iio + 7200 “ 7iio “ 72oo) 

2 T 


+2xr7f„„ - + o(*). 


(5.20) 


where we have introduced the one-loop difference 


Mf^ ^ Mr"^{g,g) - MT^{g,g). 


(5.21) 


Note that the squared matrix elements and Aif are of 0{e) and needs to be 

known up to 0{e‘^). Using the explicit results for the anomalous dimensions Eq. (5.20) 
translates into 


= -^CaMI^ {Ca - 2 NfTr) + 1 [ - ^(3 Ca Mf^ + 2 ClM[ 

- 5 Ca NfTrMI^ + 3 Cp NrTrM^^) + 4 Ca NrTr ln(-^)7W^^ 

^ ^Ca(-66 Adf® - 6Mf + CaAIo"(- 37 + 27r2)) 




18e 


+ 2NFTR{12Mf^ -9 CfM^o^ + 67Wf - 2 Ca2Wo'(-11 + ^ )) 
//2 1 1 

- 36CAA4f®ln(-^) + -CA{Mf^ - Mf) + C(e). 
s J 3 


(5.22) 


We have checked our prediction Eq. (5.22) with the explicit calculation of the gluon form 
factor in DRED and fdh [31] and we have obtained full agreement. This was of course to 
be expected, as we have verified in Section 4 that the extraction of 7 *^ from the form factor 
for H ^ gg is in agreement with its determination in SCET. 


6 Concluding remarks 

With the results presented in this paper we complete the understanding of the scheme 
dependence of IR divergent NNLO virtual amplitudes with massless particles. In particular, 
we have presented the generalization of this dependence to dred, where we have to consider 
amplitudes with external e-scalars and, hence, need the corresponding anomalous dimension 
7 e. Furthermore, we have presented a SCET approach to the scheme dependence and 
derived all anomalous dimensions again in this approach. In this way fdh and dred are 
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shown to be perfectly consistent IR regularization schemes (at least) up to NNLO, as long 
as the UV renormalization is done consistently. Concretely, this means that the various 
couplings Og, Oe and have to be distinguished. This is also the case in fdh, where at 
NNLO the only concrete modihcation appears due to the UV renormalization of the NLO 
virtual amplitudes. Our results and definitions of fdh are perfectly consistent with the 
results and dehnitions proposed in [11, 17]. 

Obviously, the virtual amplitudes are not the only ingredients needed for a calcula¬ 
tion of a physical quantity. At NNLO, also double-real and real-virtual corrections are 
to be considered. Furthermore, if there are initial state hadrons, a counterterm for the 
initial-state collinear singularities is required. All these additional contributions are also 
regularization-scheme dedendent and only once all parts are combined to a physical cross 
section, the regularization-scheme dependence cancels. 

In virtually all NNLO calculations of cross sections completed so far, CDR has been 
used. The results presented in this paper allow for using any of the other regularization 
schemes for the calculation of the virtual corrections. Using a scheme different from CDR 
often facilitates the use of efficient calculational techniques for loop amplitudes. The results 
can then be translated to obtain the virtual corrections in cdr and can be combined with 
the additional parts mentioned above, obtained again in CDR. 

Of course, it is not imperative to treat the additional contributions (i.e. the contri¬ 
butions other than the NNLO virtual corrections) in cdr. Also for these terms other 
schemes might offer advantages. In fact, a modification of a subtraction scheme at NNLO 
to the HV scheme has been presented recently [63], resulting in a reduction of the algebraic 
complexity. 

The question of the scheme (in)dependence of a full cross section at NNLO becomes 
particularly transparent if the calculation is performed in a SCET inspired way. Following 
ideas of the slicing method [64] and the g-r-subtraction method [65] , the cross section is split 
into two regions, a ’hard’ region and a ’soft’ region. In the hard region not all radiation 
in addition to the final state under consideration is soft (or collinear). At least one of the 
emitted gluons is hard. Here we are effectively dealing with a NLO calculation of a process 
for a hnal state with an additional parton and the scheme independence of cross sections at 
NLO is well established [4]. In the soft region all additional radiation is soft (or collinear) 
and a true NNLO calculation is required. For this part a SCET approach is used. This 
idea has first been applied to the decay of a top quark [32] t ^ W bX where the invariant 
mass of the jet b + X has been used for the split. Recently, the N-jettiness event-shape 
variable has been used to obtain a similar setup for differential NNLO calculations of Higgs 
plus jet [33], W plus jet [34] and Drell-Yan production [35]. 

In the soft region, the cross section factorizes into a product of hard-, soft- and jet 
functions (and beam functions if there are initial-state hadrons). The corresponding bare 
functions are all IR divergent and scheme dependent. However, we have shown that the 
properly IR subtracted soft function Sfin, Eq. (3.27), and jet functions and 
Eqs. (3.39) and (3.45), are not only finite but also scheme independent, at least up to 
NNLO. The same holds true for the hard function [66, 67] that is closely related to Adgn, 
Eq. (2.5). Hence the cross section in the soft limit can be expressed in terms of these IR 
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subtracted quantities in a manifestly scheme-independent way. 

The soft function that is required for the processes mentioned above is not the soft 
function for Drell-Yan or Higgs production as we have computed. However, the procedure 
to perform the IR subtraction (or UV renormalization in SCET language) consistent with 
the regularization scheme nsed in the computation of the bare soft fnnction is exactly the 
same. 

Since the soft, hard and jet functions are separately scheme independent, it is possible 
to nse different schemes in the computation of the various parts contributing to the cross 
section. For example, the calculation of the virtual corrections (i.e. the hard function) in 
FDH, where the helicity and unitarity methods are applicable, can easily be combined with 
the soft or jet function computed in CDR. We are convinced that this flexibility will be 
very benehcial for farther developments of fully differential NNLO calculations. 
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A Explicit expressions for the soft and jet functions 


In this appendix we give the explicit results for several quantities as a perturbative expan¬ 
sion. We use the conventions specified in Section 3.1. For most resnlts it will be snfficient 
to expand a quantity X in as and ae and write, instead of Eq. (2.1), 




as 

Att 


/CV 

Vdvr 


■ Yi . 


(A.l) 


As in Eq. (3.5) we will use the short-hand notation Xmn = and Xmn = 

The explicit resnlts for scheme-dependent quantities will be given in the 
fdh/dred scheme but we can obtain the corresponding coefficients in the hv/cdr scheme 
as Xfnn — lim^v^ _^o X^^ • 


A.l Soft functions 


It is convenient to solve the RGEs for the soft functions in Eq. (3.24) order by order in a*. 
By using the expansion coefficients of the anomalous dimensions in Eq. (3.8) one obtains 
the following scheme independent result 
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and the one and two-loop non-logarithmic coefficients have the expressions 
^2 


cY = Cr 


TT 


C 2 — <~'R 




2 vr^ 


81 


54 


TT 


+ 


4(3 _ ^ 

9 27 


(A. 2 ) 

(A.3a) 

(A.3b) 


(A.4a) 


328 \ 

w) 


(A.4b) 


The result in Eq. (A.2) is in agreement with previous calculations in [43, 47]. 

A.2 Quark jet function 

Here we list the explicit two-loop coefficients entering Eq. (3.32): 
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After renormalization and setting e —)■ 0 we obtain a finite and scheme independent quark- 
jet function. The terms containing cancel and we are left with only ols dependent terms. 
In Laplace space the quark-jet function reads 
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and is in agreement with previous results [28]. 

A.3 Gluon jet function 

Here we list the explicit two-loop coefficients entering Eq. (3.43): 
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After renormalization and setting e —)■ 0 we obtain a finite and scheme independent gluon 
jet function. The structure in Laplace space is the same as for the quark jet function, 
Eq. (A.6), 
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A.4 e- scalar jet function 


The results in this subsection depend on 04 ^ as well as Ug and Oe- We start by listing the 
explicit two-loop coefficients entering Eq. (3.48). 
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The expression for the renormalized e-scalar jet function in Laplace space is considerably 
more complicated than the corresponding expression for the quark- or gluon-jet function. 
Contrary to the quark- and gluon-jet function, there is still a dependence on Ue and a 4 e. 
The finite e-scalar jet function is given by 
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where we have kept all terms of 0{a^g), 0{a^), O(aj^) and 0(asae}, that appear in the 
structure of the equation, even if they are zero. The limit —>■ 0 has been taken and as 
usual we indicate this in the notation by dropping the bar, e.g. /3'^ = limjY^^o The 
coefficients of the anomalous dimension of the e-scalar jet can be read off Eq. (3.51). In 
particular = 0. The coefficients of the cusp anomalous dimensions can be read off 
Eq. (B.ld) and only Tioo and r2oo are non-vanishing. 

The non-logarithmic terms of Eq. (A. 12) read 
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B Anomalous dimensions 


In this appendix we collect all results for the anomalous dimensions relevant for this work 
without distinguishing the various 04 ^, 1 - 

We give the explicit results with Tr = 1/2 in the fdh/dred scheme, see Eqs. (2.20) 
and (2.21) for definitions and relations. The cdr/hv results are obtained by setting Ag = 0. 
Of course, 7 ^ is only meaningful for DRED. 
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where 0{a^) stands for a generic coupling a G {os, o;,, a 4 e,i}. 
For the /3 functions we have 
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A more complete list of coefficients for the /3 functions can be found in Ref. [17]. 
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